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ABSTRACT 



Context. Existing SWAS observations and future HIFI/Herschel data require a clear sense of the information content of water emission 
and absorption lines. 

Aims. Investigate wether the ground-state transition of ortho-IUO (lio— » loi) at 557GHz can be used to measure the column density 
throughout an interstellar cloud. 

Methods. We make use of a multi-zone escape probability code suitable for the treatment of molecular line emission. 
Results. For low abundances, i.e., X(H 2 0) < 10~ 9 , the intensity of the 1 10 — »loi transition scales with the total column of H 2 . However, 
this relationship breaks down with increasing abundance, i.e., optical depth, due to line trapping and, - for T(j ust >25K, X(H2O)<10~ 8 
and n~10 4 crrr 3 , - absorption of the dust continuum. 

Conclusions. An observed decline in intensity per column, expected if H 2 is a surface tracer, does not necessarily mean that the 
water is absent in the gas phase at large columns, but can be caused by line trapping and subsequent collisional de-excitation. To 
determine the amount of water vapour in the interstellar medium, multiple line measurements of optically thin transitions are needed 
to disentangle radiative transfer and local excitation effects. 

Key words. ISM: molecules - radiative transfer 



1. Introduction 

The lau nch of the Submill imeter Wave Astronomy Satellite 
(SWAS;|Miinick et al.l l2000) made it possible to observe emis- 
sion from the ground-state transition (lio— >loi) of ortho-H^O 
and its isotopomer ortho-H^O, to determine the abundance of 
water vapour, i.e., column density N _h,o, in various regions 
in the interstellar medium, e.g., dense and diffuse interstel- 
lar gas clouds, circ umstellar envelopes, p l anetary atmosphere s, 
and comets (e.g.. IMelnick et all |2000t fenell et all l2000allb[; 
lEncrenaz et aL[[1995t iBergin et all 120001 iGurwell et alj|200ol: 
IMelnick et al.l 120001) . In the future, the Heterodyne Instrument 
for the Far Infrared (HIFI) on Herschel will observe even 
more transitions of ortho- and para-H20 (0/P-H2O) like 2 12 — >loi 
(1669.904GHz), 2 2i ^2 12 (1661.015 GHz), 3 3-»2i 2 (1716.774 
GHz), 3i2^3(n (1097.357 GHz), 3 2 i->3i 2 (1162.910 GHz), 
ln-»Ooo (1113.342 GHz), 2 2->ln (967.924 GHz), 2„^2 02 
(752.029 GHz) and 2 20 ^2„ (1228.801 GHz); some in absorp- 
tion while others in emission. 

One of the main goals of the SWAS mission is to determine 
where and whether H 2 and 2 are the major reservoirs of oxy- 
gen through the interstellar medium. SWAS observations have 
determined the gaseous water abundance in warm dense gas 
(T > 300 K and n(H 2 ) > 10 3 cm 3 ) to be 10~ 5 relative to H 2 , 
in good agreement with chemical models for such conditions. 
However, in cold (T < 30 K), dense clouds the abundance of 
gaseous water is ~100 to 1000 times below the predictions of 
cold-cloud gas-phase chemical models. It has been suggested 
that - toward cold clouds - gaseous H 2 exists only near the 
cloud surface. Indeed, closer to the surface than an Ay of a 
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few mag, H 2 is photo-dissociated by the ambient galactic UV 
field. Deeper into the cloud, i.e., Ay of 4 ~8 mag (depending 
on density and UV intensity), H 2 Q may rapidly deplete onto 
dust g rains dBergin et alJl2000bt Icharnley et alJl200lt IViti et all 
2001). Although the derivation of the column density from ab- 
sorption observations is straightforward (column de nsity is sim- 
ply pr oportional to the optical depth in the line, see lPlume et al.l 
d2004l) ) this is not the case for emission observations. The analy- 
sis to determine the H 2 abundance now crucially depends on 
the physical properties of the gas through the collisional rate 
coefficients. Therefore, accurate constraints on the gas densities 
and temperatures are needed. 

Extensive SWAS observations of the Orion A molecular cloud 
show that gaseous H 2 correlates with C N, a surface tracer, 
rather than with C ls O, a volume tracer dMelnick & Bergin 
120051) . This result has been interpreted as evidence that gaseous 
water resides only near the surface. However, caution is needed 
when relying purely on single transition observations to draw 
such a conclusion in view of the complex rotational level struc- 
ture of the H 2 molecule. In particular, there is a fundamental 
difference between an optically thin and an effectively optically 
thin line. The latter case implies a strong coupling between line 
photons and water molecules. A full radiative transfer calcula- 
tion is needed to address this problem, since the observed inten- 
sities of molecular emission depend on a complex competition 
between radiative and collisional processes. Moreover, the ex- 
citation of H 2 also differs from that of other molecules, since 
both collisions and in frared radiation from warm dust influence 
the level populations (Takaha shi et al.|[T983l) . 
The intent of this work is to show that it is not straightforward to 
retrieve accurate information, e.g., column density, from single 
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transition observations of H2O due to the complex level structure 
of this molecule. 



2. Basic model description 

The results present ed here were obtained by application of the 
nu merical code of IPoelman & SpaanH d2005l) . described further 
in lPoelman & Spaans d2006l) . The interested reader is referred 
to these papers for a description of the underlying algorithms. 
The radiative transfer of 0/P-H2O is solved by means of a multi- 
zone escape probability method in three dimensions. By using 
a multi-zone formalism, the medium is divided into different 
zones, i.e., gridcells, each with a value for the abundance of the 
species (e.g., H2O), the density of the medium, and the temper- 
ature of gas and dust. Besides this, the cloud is characterized by 
a fixed total column density. The statistical equilibrium equation 
for a multilevel sytem can be written as 



Table 1. Model parameters 



'YjjuAij + nriijij + riiBijJ) = ^(«H,y,y + niB u J), (1) 



»j 



KJ 



where the radiative de-excitation is given by the Einstein A coef- 
ficients, By are the Einstein coefficients for absorption and stim- 
ulated emission, y (J are the collisional rate coefficients, 7 the 
mean photon radiation field, «, the population density of the zth 
level, and n the density of the medium. The level populations de- 
pend on the radiation field while the radiation field depends on 
the level populations everywher e. By introducing the concept of 
a multi-zone escape probability (IPoelman & Spaansl2005l 12006) 
the statistical equilibrium equations are decoupled from the ra- 
diative transfer equations. Eq.[TJcan be rewritten as 

^JniAijBiTij) + nnflij) = 2^(n«,-7y), (2) 



»j 



because the net absorptions, corrected for stimulated emission, 
are equal to those photons that do not escape, i.e., 



(njBji - //;#,;),. ,7 = - B(Tij)]A u , 



(3) 



i>j 



where B(t) is the probability that a photon formed at optical 
depth t in a certain direction escapes the cloud along that di- 
rection. Therefore, 



(4) 



Note that because of the large Einstein A coefficients of the wa- 
ter molecule critical densities are of the order of 10 8 -10 9 cm 
in the optically thin case. 

The background radiation field P(v, ; ) in the code consists of two 
terms: the 2.7K microwave background and the infrared emis- 
sion of dust at a temperature 7^. This is, 

P(vij) = B{v ih T = 2.7 K) + (1 - e- T ^)B{v u , T d ) (5) 

The intensity of transition i— >j (i>j) is then given by 



'//total 
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with 



Model 


n(H 2 ) 


size 


X(H 2 0)" 


N(H 2 ) 6 


T„as 


Tdust 




[cm" 3 ] 


[pc] 




[cm" 2 ] 


[K] 


[K] 


I 


10 4 -10 6 


0.002-0.2 


io-w-io- 6 


6xl0 21 


50 


no 


II 


10 4 -10 6 


0.002-0.2 


10- 10 -10- 6 


6xl0 21 


50 


15 


III 


10 4 -10 6 


0.002-0.2 


io- lo -io- 6 


6xl0 21 


50 


25 


IV 


10 4 -10 6 


0.002-0.2 


io- lo -io- 6 


6xl0 21 


50 


50 


V 


10 4 -10 6 


0.002-0.2 


io- lo -io- 6 


6xl0 21 


30 


15 



" Abundance of H 2 0; Total column density through the centre 

S(vy) is the source function at frequency Vy. 
Collisional rat e coefficients fo r inelastic collisions between o/p- 
H2O and He (iGreen et alJll993 j ), and for collision s between o/p- 
H 2 and both o-H 2 and p-H 7 ( iPhillips et alj[l996l) are adopted. 
We adopt the expression for the ortho-to-para ratio (OPR) of H2, 
in thermal equilibrium, defined by 



OPR = 



(2/ + 1) 2(27 + l)exp(- £ ° (y ^ ) ) 
(27 p + 1)2(27 + l)exp(- Ep( ^ c) ) ' 



(8) 



Ayiocal = niAijhVijB(Tij){[S (Vij) - P(Vij)] /S (vy)}, (7) 



where 7 D and 7 p are the total nuclear spin, corresponding to 
wether the hydrogen nuclear spins are parallel (7 Q = 1, Tf) or 
anti-parallel (7 p = 0, TD- The sum in the numerat or (denomina- 
tor) ex tends over all ortho (para) levels (7, K a , ATA lMumma et al.l 
( 1987). The code has been tested extensively against (analytical) 
benchmark problems presented at the radiative transfer w ork- 
shop held in Leiden (2004), see IPoelman & Spaans! d2006). It is 
found that the level populations are completely consistent with 
th e solutions of other Mon te Carlo and ALI codes, as presented 
in lvan der Tak et all d2005l) 



3. Model results 

We calculate, in the case of a homogeneous sphere and as a func- 
tion of impact parameter, the surface brightness for the 1 10— > loi 
ground-state transition of ortho-H^O. The density and abun- 
dance of H2O are the main parameters of the model. All the 
models have a constant total column density NCH2) of 6 x 10 21 
cm -2 through the center of the cloud, corresponding to a total 
Ay of ~3 mag. As a result, the physical size of the cloud is 
inversely proportional to the density of the medium. The den- 
sity ranges from 10 4 to 10 6 crrT 3 covering the relevant range for 
dense molecular clouds such as the Orion ridge. The tempera- 
tures, ranging from 30 to 50K for the gas and from 15 to 50K for 
the dust, were chosen to represent the mean observed tempera- 
tures towards star forming molecular clouds such as the Orion 
ridge which have been the focus of the SWAS effort. In model I 
we ignore the emission from dust and CMB, whereas in model 
IV a dust temperature of 50K is assumed, in order to assess the 
influence of the dust and temperature, respectively, on the exci- 
tation of the water molecule. The temperatures in model V are a 
factor of ~2 lower than in the other models since to maintain a 
temperature of the gas of 50K throughout the cloud one needs a 
strong UV radiation field, which is not always the case. Note that 
the gas and dust temperatures are independent of cloud depth. 
The intent of the paper is to illustrate the excitation and radiative 
transfer effects assuming a 'simple' cloud model, not to model 
a realistic cloud. A Galactic dust-to-gas ratio of 10~ 2 by mass is 
assum ed. We adopt the dust opacities of lOssenkopf & He nning 
(1994) (Column 5 of their Table 1). In all the models, except 
for model I, the dust optical depth Td us t through the center of the 
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Fig. 1. The intensity of the ortho-H^O ground state transi- 
tion for a homogeneous sphere with H2 densities of I0 4 (top), 
I0 5 (middle), I0 6 (bottom) cirT 3 and temperature of the gas of 50 
K. The dust emission as well as the CMB radiation are ignored 
{i.e., model I). In every case, the total column density (H2) is kept 
constant. Lines are plotted for an abundance of H2O, relative to 
H2, from 1CT 10 (upper line/red) to 1CT 6 (lower line/light blue) as 
function of H2 column density along the line of sight, where 2 X 
10 19 cm~ 2 is at the edge, and 6 x 10 21 crrr 2 through the center 
of the cloud. F-axis in units of erg s sr 1 . The position of the 
t=1 surface is displayed with a cross for X(H2O)=l(T 9 -10~ 6 . 



Fig. 3. The intensity of the ortho-H20 ground state transition in 
case of a homogeneous sphere with densities (H2) of I0 4 (top), 
\0 5 (middle), I0 6 (bottom) cirT 3 and temperatures for gas and 
dust of 30 K and 15 K, respectively (i.e., model V). In every 
case, the total column density (H2) is kept constant. Lines are 
plotted for an abundance of H2O, relative to H2, from 10~ 10 (up- 
per line/red) to 10~ 6 (lower line/light blue) as function of H2 col- 
umn density along the line of sight, where 2 x 10 19 crrr 2 is at the 
edge, and 6 x 10 21 cirT 2 through the center of the cloud. F-axis 
in units of erg s sr~' . 
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1 x1 2 2x1 2 3x1 2 4x1 2 5x10 2 6x10 2 
COLUMN DENSITY (H 2 ) 

Fig. 2. The intensity of the ortho-H20 ground state transition for 
a homogeneous sphere with density (H2) of 10 4 cm -3 for model 
II (top), III (middle) and IV (bottom). Lines are plotted for an 
abundance of H2O, relative to H2, from 10~ 10 (upper line/red) 
to 10^ 6 (lower line/light blue) as function of H2 column density 
along the line of sight, where 2 x 10 19 crrr 2 is at the edge, and 
6 x 10 21 cm -2 through the center of the cloud. F-axis in units of 
erg s _1 sr _1 . 



cloud at the frequency of the ground state transition of 0-H2O, 
i.e., t v= 556.936ghz or T / i = 538 jUm , is 10~ 3 . Within each model the wa- 
ter abundance, X(H20), ranges from 10~ 10 to 10~ 6 . The param- 
eters for the different models are shown in Table Q] Throughout 
the models, a velocity dispersion of 1 km s~' is adopted, typical 
for a cloud with moderate turbulence. Note that the results pre- 
sented in this paper depend on the adopted velocity dispersion. 



0.0001 



0.0010 



uj 0.0010 



10' 6 



10' 7 



RADIUS (cm) 



Fig. 4. Shown here is the population of the lot level of 0-H2O as 
function of radius. Top, middle, bottom panel are the results in 
case X(H 2 0) = 10" 10 , 10~ 9 , 10~ 8 , respectively. Dashed, dotted 
and dashed-dotted lines represent the results of model I, III and 
IV, respectively. 



A higher (lower) velocity dispersion will decrease (increase) the 
optical depth for a given transition, hence having an impact on 
the excitation of the molecule. 

Fig.[T]-[3]present the basic results of this work. We plot, as a func- 
tion of N(H2), i.e., different impact parameter, the 1\q— »1oi line 
intensity above the continuum per unit column density of 0-H2O, 
see Eq.|6] The results for model I, in which we ignore the contri- 
bution of dust and CMB on the excitation of the water molecule 
are displayed in Fig. Q] The position of the t-1 surface is de- 
noted by a cross for X(H 2 O)=10- 9 -10- 6 . Models II, III and IV 
behave in a similar way as in model I in the case n(H2)>10 5 cm~ 3 . 
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10"'° 10~ 9 10~ 8 10~ 7 10~ 6 

ABUNDANCE H 2 



Fig. 5. Conversion of the results of model I {top) and IV {bot- 
tom) into surface area weighted values. In each panel the dashed- 
dotted, dotted, and dashed curve represent the results in case the 
density n(H2) is 10 6 , 10 5 and 10 4 cirT 3 , respectively. F-axis in 
units of erg s _1 sr _1 . 

For this reason, we only plot the outcome in the low density case 
for models II, III and IV in Fig. [2] The results for model V are 
plotted in Fig. [3] 

The following trends can be identified, which will be discussed 
in Sect.|4] 

First, in the case of n(H 2 ) > 10 5 cm 3 and X(H 2 0) < 1(T 8 , i.e., 
t<10, a linear relationship holds between the number of pho- 
tons escaping the cloud and the impact parameter, i.e., I/Nh 2 o is 
constant. However, this relationship breaks down at high opti- 
cal depth, i.e., H2O > 10~ 8 , for all the models. Second, in all 
the models, except in model I, absorption occurs in the low 
density case when abundances are low, i.e., X(H20) < 1CT 8 , 
which becomes more apparent as the dust temperature increases. 
However, the amount of absorption is moderate as the self- 
reversal in the center of the line is small. Third, when the H2O 
abundance exceeds 10~ 7 , in all the models, the ratio of the in- 
tensity to the column density decreases near the edge of the 
cloud (N(H 2 )~5xl0 20 ctrT 2 ). Fourth, for high optical depth, i.e., 
X(H 2 0)>1(T 7 , and n(H 2 )=10 6 cm" 3 , I/N H2 o decreases with in- 
creasing column density. Fifth, lowering the gas and dust tem- 
peratures by a factor of ~2 (model V) does not lead to signifi- 
cant differences in the shapes of the curves. That is, models II-IV 
and model V experience similar complicating radiative transfer 
effects. 

4. Discussion 

The asymmetry of the water molecule causes the rotational lev- 
els to split into a number of different ladders, so called 'K- 
ladders', characterized by different values of the projection of 
the angular momentum onto the principal axes of the molecule 
{Jk_k + )- Radiative transitions occur rapidly between levels in 
each ladder but are much slower between levels in different lad- 
ders. This leads to a spectrum more complex than linear or sym- 
metric top molecules, e.g., CO, NH3. Hence, it is not straight- 
forward to disentangle the different processes that contribute to 
the observed spectrum. We now describe the different effects that 
play a role in the interpretation of the figures. In this, the 'edge' 
and 'centre' of the cloud refers to an impact parameter of 1 and 
0, respectively. Note that the use of spherical models leads to an- 



gular non-trivial re-distribution of line photons. The same holds 
for continuum photons within the line profile frequency range. 
First, one can see in Figs. Q] and [3]that the curves as function of 
total column density (H 2 ) are constant for X(H 2 0) < 10~ 8 and 
n(H 2 )>10 5 cm~ 3 . In this limit, collisional de-excitation and scat- 
tering effects are negligible. Eventually every photon produced 
in the cloud will escape the cloud with few interactions with 
the surrounding medium. Note that the number of scatterings N 
to escape depends on the optical depth. In this regime r <k 1, 
therefore few photons are scattered and N ~ t. With increasing 
optical depth, i.e., X(H 2 0) > 10~ 8 , more effects have to be taken 
into account. All models show a drop near the edge of the cloud. 
Because of the increasing optical depth, line-scattering effects 
become important. Thus, line photons then tend to escape in the 
direction with the lowest optical depth rather than tangentially 
to the cloud surface, causing the dip near the edge. However, to- 
wards the centre of the cloud, the optical depth increases with 
orders of magnitude. The photons will undergo numerous scat- 
terings for t » 1, with N ~ t 2 , and eventually will escape in the 
line wings. 

Second, at densities as low as 10 4 cirT 3 , and abundances not 
exceeding 1CT 8 (modest optical depth), the line is strongly sub- 
thermally excited, and radiatively colder than the dust back- 
ground. Hence, the line appears in absorption. The decrease in 
intensity/N(H 2 0) shown in Fig.|2](red and green curves) now in- 
dicates that lines of sight through the cloud center are no longer 
contributing evenly to the emissivity across their entire column. 
The line is not strictly in absorption yet, but it has developed 
an intensity dip around line center. Thus, the presence of dust 
causes the trend in the intensi ty per column in this regime to de- 
crease (Takahashi et al. 1983). This behaviour is not seen in case 
n(H 2 ) > 10 5 cirT 3 , as in this regime collisions are the dominant 
process in the excitation of the water molecule, thereby nulli- 
fying the effect of dust emission. The influence of dust on the 
excitation/level populations of water is plotted in Fig. [4] where 
the relative population of the 1 10 level is displayed. One notices 
that for warmer dust level 1 10 is more populated. In essence, dust 
continuum emission will tend to drive the level populations to- 
wards a Boltzmann distribution at the temperature of the dust. 
For a given density, the effects of radiative excitation by dust 
continuum emission is more pronounced for higher dust temper- 
atures (e.g., higher continuum intensities). 
Third, the effect of photon trapping is to lower the density at 
which LTE is approached, i.e., after each absorption, the gas has 
a chance to collisionally de-excite the species and return the ex- 
citation energy to the thermal bath of the gas, see Eq. [4] For 
optically thin gas, the critical density of the ground state tran- 
sition at 557 GHz of o-H 2 is ~10 8 cm 3 at 50 K. The optical 
depth, through the centre of the cloud, varies from 0.1 to 10 3 
when the abundance rises from 10 -10 to 10 -6 in all the mod- 
els. Hence, for high abundances, i.e., X(H 2 0) >10~ 7 , the effec- 
tive critical density drops to 10 5 -10 6 cirT 3 , since for high optical 
depth (3{t)~ 1 / t. Collisional de-excitation processes then become 
important in the regime where n(H 2 ) > 10 5 cm 3 , and X(H 2 0) 
> 10~ 7 . It is seen in Fig. Q] and [3] that for n(H 2 )=10 6 ctrT 3 and 
X(H 2 O)=10~ 6 the I/Nh,o drops as function of impact parameter. 
In this part of parameter space collisional de-excitation is impor- 
tant and the probability that line photons are lost to the thermal 
bath through collisional de-excitation during one of the many 
scattering events is high. 

Fourth, calculations are performed for model V (Fig. [3} with gas 
temperatures a factor of 2 lower relative to the temperatures used 
in model II. We find that the shape of the curves are not affected 
by such a change in gas temperature. However, it affects the dis- 
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tribution of the level populations and thus the absolute intensity 

in the lines. Hence, temperature variations cannot dispense of 

the radiative transfer effects studied in this work. 

To summarize, we plot in Fig.[5J as a function of abundance, the 

average intensity emanating from the cloud for model I and IV, 

i.e., 

$ h m ^Jb)2nb db 
f N H2 a(b)2nb db ' 

with b the impact parameter. One notices a drop by a factor of 
~2-5 in case n(H2) is 10 5 — 10 6 crrT 3 , respectively. Note that 
with the assumption of an effectively optically thin line one 
would underestimate the water column by these same factors. 



5. Astrophysical implications 

The intensity of the ground-state transition of 0-H2O is driven by 
a combination of the ambient gas and dust temperatures on the 
one side and by the density of the surrounding medium on the 
other side. It is this interplay, together with the complex struc- 
ture of the molecule that drives the level populations. To interpret 
existing SWAS and future HIFI data a clear sense of the infor- 
mation content of the water lines is needed. 
SWAS observations of the lowest rotational transition of o-H^O 
of the Orion A molecular cloud show that gaseous water corre- 
lates much better with the near surfac e tracer CN than with th e 
volume tracer C 18 0, as presented in Melnick & Bergin (2005). 
Through these observations - in which it is assumed that the 
ground-state transition of ortho-H20 is effectively optically thin 
- one concludes that water is a surface tracer. This is plausi- 
ble from a chemical point of view in which photo-dissociation 
destroyes the water molecule near the surface. Further inwards 
the cloud the water abundance reaches its equilibrium value 
through photodesorption of F^O-ice and photodestruction of 
H20-gas until it freezes-out onto dust grains deeper into the 
cloud. However, we have shown, as seen in Fig. [5] that for 
r>10 and n>10 5 crrT 3 the effectively optically thin assumption 
no longer holds. Hence, under these conditions one is limited 
to observing the 't=10' surface and cannot use water to trace 
the cloud's volume, even when it is present, i.e., not frozen out 
(Cernicharo, private communication). Therefore, as CN is a sur- 
face tracer and the water intensity originates from a layer of gas 
with an optical thickness of 1 - depending on local excitation 
conditions this layer is a surface layer - the CN intensity corre- 
lates much better with the H2O intensity and not with the vol- 
ume tracer C ls O Thus, the anti-correlation of H2O with C ls O is 
partly due to optical depth effects, and is not neccesarily a result 
of chemical changes. As a consequence, the presence of water 
past the t=1 surface cannot be ruled out. 

We would also like to point out here that the most interesting 
aspect of the correlation of the water line intensity with the CN 
line intensity is the fact that both are observed to vary by a fac- 
tor ~100. Theoreticall y, the CN abundance is ex pected to scale 
with density squared (Boger& Sternberg 2005), indicating the 
importance of density variations over the Orion molecular cloud. 
Given the results presented in this paper, we surmise that these 
density variations will hamper the interpretation of the water ob- 
servations. 

In order to deduce the total water column along the line of sight, 
additional information is needed from other - effectively opti- 
cally thin- lines, which will be observed with future missions 
such as Herschel/HIFI. 
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